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Abstract 

O 

The expansion dynamics of hot electron-positron-photon plasma droplets is dealt with within 



Oh 



relativistic hydrodynamics. Such droplets, envisaged to be created in future experiments by ir- 
radiating thin foils with counter-propagating ultra-intense laser beams, are sources of flashes of 



gamma radiation. Warm electron-positron plasma droplets may be identified and characterized by 
^ a broadened 511 keV line. 
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I. INTRODUCTION 



The rapid progress in ultra-intense laser technology pQ triggered the interest in properties 
of relativistic plasmas consisting mainly of electrons (e~), positrons e + and photons (7). For 
instance, in pHH] various processes, including thermal (kinetic) and chemical equilibration, 
have been considered for a hot charge-symmetric electron-positron-photon (e~e + 7) plasma, 
while the kinetics of a temporarily created electron-positron plasma has been analyzed in 
[6]. Hot plasmas, may be charge-symmetric or asymmetric or with ion load, have inter- 
esting features [7], e.g. with respect to collective excitation modes (including the van Hove 
singularity), astrophysical situations and secondary particle (such as muon or pion) produc- 
tion. Furthermore, an optically thick e~e + 7 plasma represents the Abelian analog of the 
non-Abelian quark-gluon plasma [8] currently investigated at large-scale accelerators. 

The scenario for droplet creation can be visualized in the following way: Suppose that 
the radiation pressure of counter-propagating ultra-intense laser beams on a thin foil is 
so high that a substantial compression accompanied by sizeable pair creation is achieved. 
Estimates for such a process have been presented in [9]. Extrapolating these estimates 
one may envisage the idealized scenario of a hot and optically thick e~e + 7 plasma, with 
temperatures exceeding the 1 MeV scale. Optically thick e~e + 7 plasmas play an important 
role in the gamma ray bursts phenomena [10J. 

In [3] it has been pointed out that an e~e + 7 plasma droplet at initial temperatures T 
much larger than the electron mass m emits an intense gamma ray burst. However, the 
expansion dynamics was considered in [3] in a schematic manner. Nevertheless, the time 
structure and the spectrum of the burst depend crucially on the expansion dynamics of the 
plasma droplet driven by the thermodynamic pressure gradients. It is necessary to consider 
a realistic expansion dynamics as well as the particle production rates in the relativistic 
framework for a proper understanding of e~e + 7 plasmas produced in the laboratory and 
astrophysical sites such as gravitational collapse. The aim of the present paper is to quantify 
the droplet expansion dynamics by means of the relativistic hydrodynamics. This provides a 
realistic evaluation of the time evolution of the temperature profile and velocity profile. The 
latter one is important for properly accounting for the blue/red shift effects of the emitted 
photons. 

Another avenue towards the creation of e~e + plasma droplets has been envisaged in [B] 
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as a result of the dynamical Schwinger process. Due to the short (sub cycle) time scales the 
thermalization and the collective dynamics are not an issue. Instead, the emergence of e~e + 
annihilation into 7 pairs may be employed as diagnostic tool. 

Our paper is organized as follows. In section |TT] we consider the spherically symmetric 
expansion of plasma droplets within the relativistic hydrodynamics. This input is used in 



section III to estimate the spectrum of emitted photons. Section IV gives an account of 
the spectral emissivity for a warm plasma with a temperature in the order of the electron 
mass and focusses on the onset of the broadening of the 511 keV line when achieving higher 
temperatures. Section [V] is devoted to the summary. Appendix [A] presents the expansion 
dynamics with planar geometry. 



II. RELATIVISTIC HYDRODYNAMICS 

A sufficiently large e~e + 7 plasma droplet can be dealt with by means of relativistic 
hydrodynamics. The equations of motion imply local energy and momentum conservation 
expressed as 

= 0, (1) 

where fi, v — 0, . . . , 3 denote Lorentz indices and a colon stands for the covariant derivative. 
For an isotropic medium (gas or fluid or plasma) the energy-momentum tensor reads 

= (e + p)u"u v - pg^ + ■■■ . (2) 

Here, the four-velocity is u M with normalization u^u^ = 1 and g^ v denotes the metric ten- 
sor, diag(g Ml/ ) = (1,-1,— 1,-1) in Cartesian coordinates. The energy density is e, and p 
stands for the thermodynamic pressure. The ellipsis • • • in (|2j) denote dissipative terms 
being proportional to transport coefficients, such as shear and bulk viscosities, for instance. 
Neglecting the latter ones, ([TJ [2]) can be solved once an equation of state is given. For 
T > m, the equation of state of a relativistic ideal gas can be employed, that is e = 2>p. 
We consider a charge-symmetric plasma with zero net density. Accordingly, the chemical 
potential is zero, too. The pressure and temperature for our equation of state are related via 
p = deff^T 4 (we employ units with h — c—1) with d e ff = 2 + |4 for photons (first term) 
and electrons and positrons in equilibrium. This is the leading-order term in an expansion 
in powers of m/T. 
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Special symmetries simplify Q tremendously. We consider here spherically symmetric 
droplets (planar slab geometry is considered in appendix A). Then (JIJ degenerates to two 
coupled equations, one for the pressure gradient and one for the radial velocity gradient. 
These equations can be solved by the methods of characteristics (cf. [H] for details). Nu- 
merical results are exhibited in Fig. 1. 
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FIG. 1: Temperature (left) and velocity (right) profiles at various time instants for the spherically 
symmetric expansion. 



The solution for the given equation of state can be represented for scaled quantities, T/T 
and r/R, where To and R stand for the initial temperature and radius, respectively, of the 
e~e + 7 droplet. We assume initially homogenous droplets, i.e. T(t = 0,r) = T for r < R 
and T(t = 0, r) = for r > R. The initial radial profile is v (t = 0, r) = for r < R. 
Initially, a rarefaction wave travels towards the center with the velocity of sound v s = c/y3 
(corresponding to a Riemann type solution) and becomes reflected then, causing the non- 
trivial radial temperature profile with somewhat lower temperature in the center. We allow 
the temperature evolution until (T) ~ O.lTo. For late times, the temperature drops further, 
but our assumed equation of state as well the photon rate below based on the perturbative 
estimates do not longer apply. These late times until disassembly are not important for 
the hard photon emission off the expanding droplet. The scaled time to cool from T to 
(T) ~ O.lTo is about ct = 3.5R. That means, for instance, a droplet of initial radius R = 2 
nm and To = 10 MeV with energy content 13 kJ expands within the relaxation time for 
pairs and photon chemical equilibration, 0.025 fs, and cools to (T) ~ 1 MeV. 

The radial velocity profiles in Fig. 1 expose the rapid expansion of the cool surface layers 
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nearly at velocity of light. In particular, due to the rapid expansion, the droplet cools 
faster that what was assumed in the schematic model in |3j. This is reflected in the photon 
emission considered in the next section. 



III. PHOTON EMISSION 

Having at our disposal the realistic temperature and velocity profiles, we can now proceed 
to calculate the photon emission. It should be emphasized that the leakage of energy via 
photon emission is not properly accounted for in ([Tj [2J (cf. [12J for hydrodynamic schemes 
with energy leakage). It turns out however that the emitted energy during the considered 
evolution is small in comparison with the internal and kinetic energies so that a consistent 
hydrodynamic treatment with radiation transport analog to stellar core collapse simulations 
is not needed. 

We employ the emission rate from thermo-field theory [13] (cf. figure 1 in [3]) 



to- 



d*xd 3 k (2tt; 
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(3) 



where f+i(oo) = (exp(u/T) + denotes the Fermi distribution function, and the dynam- 
ically generated asymptotic mass of the electron is = 2m 2 h . The thermal mass squared 
is determined by m\ h = e 2 T 2 /8 with e 2 = Ana and the C coefficients are given explicitly 
in [3]. = u^k^g^ takes properly into account the red/blue shifts via the scalar product 
of the medium's four-velocity u^(t,x) and the photon's four-momentum k v = (oo,k) The 
rate includes 2 <H- 2 Compton and annihilation processes (first three terms), bremsstrahlung 
processes (C&) and off-shell pair annihilations processes (C a ). 

The space-time integrated emission spectrum includes only such photons with a mean 
free path A 7 being larger larger than the distance to the cool droplet surface in a given 
direction. The mean free path is locally determined by A 7 = T" 1 with the damping rate 
T 7 (u;) = @EL e w l T d i^ dik which depends in turn on the photon energy. 

The resulting energy-weighted photon spectrum is exhibited in Fig. [2] Similar to the 
schematic model [3] it extends beyond 50 MeV with notable strength. In contrast to [3], 
however, the maximum of the spectrum is an order of magnitude smaller. The origin of this 
difference may be traced back to the improved description of the droplet dynamics which 



makes the droplet cooling considerably faster. The hard tail of the photon spectrum is 
emitted at early times and does not change at later times. In contrast, the softer part with 
ui ~ 10 MeV gets significant contributions at later times. Comparing the spectra at t — 4 
nm/c and t = 7 nm/c, one concludes that the very late radiation contribution for T < 1 
MeV will contribute to photon energies u < 5 MeV. 
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FIG. 2: Energy- weighted spectrum of emitted photons as a function of the energy for various 
times for T > 1 MeV. 



IV. BROADENING OF THE 511 KEV LINE 

The photon emission rate ^ applies for T > m, i.e. a hot plasma. In a warm plasma 
with T ~ m or T < m, one expects a significant contribution from the annihilation process 
e~e + — > 77. The corresponding inclusive photon rate for this process is obtained by the 
kinetic theory estimate 

dR M f°° ds rtmax 



Emax 



x / iE ME)M»-E) (4) 
JE min a/ (E max — E)(E — E min ) 
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with Af = A and 




(5) 
(6) 

, -B ± VB 2 - AAC 

£jmm — — , [I) 



m a ,r 



A = -s 2 , B = 2s(us + 2ujt' -vt'), (8) 
C = st'{s + f) - u 2 s{s - Am 2 ) 

- 2uust' - vH' 2 - Auism 2 s + s 2 m 2 (9) 

along with the additional kinematical restriction — As 2 t'(s + t') (1 + sm 2 /((s + t')t')) > 0. 
Ai(s,t) is the lowest-order invariant matrix element for the process e~e + —> 77 (cf. [1]), 



m 2 — u m 2 — t , m 2 m 2 



2 2 1 2 N 

m m 

+ 



m' 2 — t ' m 2 — u 



|A4( S ,t)| 2 = 327r 2 a 2 I ^^ + ^^1+4^— + 4^— -4 

m 2 — s m 2 — u m 2 — t m 2 — t 

"(10) 

with u = 2m 2 — s — t, and the distribution functions are fs{y) = (exp{y/T} + «S) _1 . For 
T ^> m, corresponding to m — > 0, one recovers the form presented in [H]; the rate scales 
oc T 2 . In the opposite region oj ^> T, where the Boltzmann approximation for the f + \ terms 
is justified [Hj, one finds 

dR M T _ W/T 

U d±x(Pk 2(27r) 5 w e 

f°° rla 

x / -( s -4m 2 ) 2 a( S ) lna-e- 8 ^))- 1 , (11) 

J Am 2 S 

where the total cross section a(s) follows from integrating da/dt = \Ai\ 2 /(16ir(s — Am 2 ) 2 ) 



within t m in — tmax- Due to the interplay of the factors e~ W//T and ln(l — e _S// ( 4a;T ^) _1 in (11) 
the rate has a sharp maximum at 00 = m for T ^ m. This is the 511 keV annihilation 
line, at very low temperatures, which broadens rapidly with increasing temperature. The 
high-frequency tail, to ^> max(T,m), behaves as codR/d 3 k oc e _tJ//T . The full width at half 
maximum is Aw = ^ —AmT In |, for T<m, with maximum oc e~ 2m//T . 

The rate Q is displayed in Fig. [3] for a sequence of temperatures. With increasing 
temperature the relative kinetic energy of the annihilating electron-positron pairs increases 
and, consequently, the photon energy increases in the medium's rest frame and becomes 
larger than m. This broadening of the 511 keV annihilation line, once isolated from other 



processes with continuum radiation such as the Compton and the various bremsstrahlungs, 
may serve as a measure of the achieved temperature in a warm e~e + 7 plasma. At T > 0.1m, 
the broadened annihilation line changes gradually into a smooth continuum with typical 
exponential tail. For T > m, the maximum of the emission rate is shifted to u> > m. The 
annihilation is then part of the complete rate ([3]). 

The presented annihilation rates apply also to the transiently created e~e + pairs by 
the dynamical Schwinger process, as envisaged in [6], when approximating the distribution 
functions by thermal ones. 
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FIG. 3: Photon emission rate from e~e + pair annihilation as a function of photon energy for 
various values of the temperature. Note the scale factors in the legend. The vertical dashed line 
indicates the position of 511 keV line. 



V. SUMMARY AND DISCUSSION 

In summary we consider the expansion dynamics of hot plasma droplets which can be 
created with next-generation ultra-intense laser beams. We present here the dynamics of 
charge-symmetric and equilibrium e~e + 7 plasma droplets by means of relativistic hydrody- 
namics. Our results support the estimate in [3] pointing to gamma flashes emitted off the 
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rapidly exploding droplets. 

Clearly, our considered scenario is still idealized. The surplus of electrons over positrons 
needs to be considered in future analysis. This can be accommodated by a finite chemical 
potential steering the asymmetry of electrons and positrons. Furthermore, at given temper- 
ature, implying kinetic equilibrium, the positrons may not be chemically equilibrated. This 
requires an additional dynamical variable, e.g. a positron fugacity supplemented by appro- 
priate rate equations. Also, ion impurities may have an impact on the plasma dynamics. 

As indicated by the largely differing expansion time scales for spherical and plane expan- 
sion patterns, the flow symmetry may play an important role. To clarify that issue, more 
simulations are needed along with realistic initial conditions. The lessons learnt from the 
investigation of the quark-gluon plasma [8], however, indicate that it is useful to develop 
parallel with experiments the theoretical description after first initializing considerations. 

Finally, we mention the broadening of the 511 keV line as useful tool for thermometry 
of small electron-positron-photon droplets produced in the laboratory and in astrophysical 
sites. 
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Appendix A: Planar Slab Geometry 

A one-dimensional planar expansion (cf . [11] ) needs a much longer time to cool the system 
down to T to (T) ~ 0.1T . This is exhibited in Figs. |l]and[5j 

To understand qualitatively this huge difference to the spherically symmetric expansion, 
one may mention that for the latter one an increase of the radius by a factor of 10 4//3 is 
needed to cool the medium from To to O.lTo when neglecting the conversion of internal 
energy into kinetic energy and emission. For planar geometry, in contrast, the thickness of 
a slab must increase by a factor 10 4 . This estimates explain the much longer time scales to 
cool the medium to a tenth of the initial temperature. Of course, the transverse expansion 
of a finite slab of matter will become important and requires more involved multi-parameter 
solutions of the hydrodynamical expansion. 
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FIG. 4: Temperature profiles at various time instants for a planar one-dimensional expansion. 
Left (right): short (long) times. 
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